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Abstract 

Exploiting final-state interactions and/or identity interference, analysis of 
anisotropic correlations of particles at low-relative velocities yields information on 
the anisotropy of emission sources in heavy-ion reactions. We show that the use of 
cartesian surface-spherical harmonics in such analysis allows for a systematic ex- 
pansion of the correlations in terms of real angular-moment coefficients dependent 
on relative momentum. The coefficients are directly related to the analogous coeffi- 
cients for emission sources. We illustrate the analysis with an example of correlations 
generated by classical Coulomb interaction. 
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Correlations of particles at low relative velocities are commonly used for as- 
sessing geometric features of emission zones in nuclear reactions with multi- 
particle final states [1,2]. At low relative momentum, the correlations exhibit 
structures, due to final-state interactions and/or identity interference, that are 
more pronounced for tighter emission zones. When particle emission zones are 
deformed, corresponding shape anisotropies are observed in the correlation 
function which can be studied as a function of the orientation of the relative 
momentum. Measured anisotropies in the correlation functions thus provide 
insight into the shape of the emission zone which is intimately connected to 
such aspects of the emission as reaction geometry, collective expansion, emis- 
sion duration and differences in emission times for different species [3,4,5]. 

Final-state interactions and identity interference link the measured correlation 
function, 7£p(q) (or alternatively Cp(q)), to the source function, <Sp(r), which 
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provides the probability that two particles of the same velocity, whose total 
momentum is P, are separated by a distance r at emission [6,7]: 



ftp(q)=C P (q)-l = 



d 6 N ab 
d 3 p a d 3 p b 



dN a dN b 
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- 1 



- 1 



Sp 6 (r). (1) 



Here, q and r are the relative momentum and relative spatial separation as 
determined by an observer in the two-particle rest frame and 0*") is the relative 
wave function for the asymptotic momentum q. For particles with intrinsic 
spins, the square I0* - - 1 ! 2 is averaged over spins, ensuring its dependence only 
on q, r and the angle between the vectors, # qr . 

Various means have been employed to analyze anisotropies in correlations 
to provide quantitative measurement of the anisotropies in a source function. 
When analyzing identical pion correlations, the most common means has been 
to present projections of the correlation function along fixed directions of the 
relative momentum while constraining the remaining momentum components. 
Such projections often neglect regions of the data outside of the employed cuts 
and they can be clumsy if one is unaware of the best directions to orient the 
coordinate system. A second means for analyzing anisotropies involves fitting 
to three-dimensional parameterizations, such as for gaussian or blast-wave 
models [8]. Since source parameterizations typically employ on the order of 10 
parameters, it can be difficult to determine the confidence of the fits due to 
the inherent complexity of cross correlations. 

A third class of approaches involves applying tesseral spherical harmonics 
Yim (regarding terminology for harmonics see [9]), to express the directional 
information in the correlation function and in the source function [10], 

n £m (q) = (47T)- 1 / 2 J dfi q r/JOJ K(q) , 

S em (r) = (47T)- 1 / 2 J dn r Y; m (Q r ) <S(r). (2) 

Here, the labels indicating the total momentum P and the particle types are 
suppressed. In terms of the projections, Eq. (1) can be re-expressed [10] as 

n tm {q) = 4nJ drr 2 K E (q,r)S Em (r) , (3) 



where 



1 f r 
(q,r) = -J d cos 6> qr P £ (cos 6> qr ) [{^(q, r, cos qr ) 
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Table 1 

Cartesian surface-spherical harmonics of rank £ < 4. Other harmonics can be found 
by permuting indices or switching indices, i.e., x <-> y, x <-> z or y <-> z. 
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^4xx = n ir ~~ 1/3 
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A x % y = n 2 x n 2 y - (l/7)n 2 x - (1/7)^ + 1/35 


-4*1, = n^n y - (l/5)7ij, 


.Aixj/* = n£ n y n z - {l/7)n y n z 



Since a given £m projection of 1Z is only related to the same £m projection 
of iS, the complexity of deconvoluting the shape information present in TZ is 
enormously simplified. 

The disadvantage of expansion in the tesseral harmonics Yi m is that the 
connection between the geometric features of the real source function «S(r) 
and the complex valued projections Se m (r) is not transparent. The Y^ m har- 
monics are convenient for analyzing quantum angular momentum, but are 
clumsy for expressing anisotropics of real-valued functions. To rectify the 
shortcoming we propose performing a similar analysis using real- valued carte- 
sian surface-spherical harmonics [11,9] that directly address moments of the 
analyzed anisotropic functions. Cartesian harmonics are based on the products 
of unit vector components, n ai n a , 2 ■ ■ -n ar Due to the normalization identity 
n 2 x + n 2 + n\ — 1, at a given £ > 2, the different component products are not 
linearly independent as functions of spherical angle; at a given £, the products 
are spanned by tesseral harmonics of rank £' < £, with £' of the same evenness 
as £. Projecting away the lower-rank tesseral components, £' < £, is a linear 
operation in the cartesian rank-£ space, represented by a rank-2£ tensor V^ 2t> 
that produces cartesian rank-£ harmonics A^\... ai (Vi): 

^■al-a e = 7 > iJ..a e ,a' 1 ...a' i n *'i " ' n <*' t ' ( 5 ) 

a' 1 ...o.' l 

Tesseral harmonics of rank £, are, conversely, expressible in terms of the powers 
of unit- vector components of the order £' < £, with £' of the same evenness as £. 
In the end, the harmonic A^ from (5) may be expressed as a combination of 
the powers of n a of order £' < £, where £' is of the same evenness as £, as shown 
in Table 1. The leading term is a simple product of unit vectors with the same 
indices as A. Since cartesian harmonics A^ are linear combinations of tesseral 
harmonics of rank £, the product r e A^ satisfies the Laplace equation. 

Because A^ is a linear combination of the tesseral harmonics of rank £, and 
because the kernel in Eq. (3) depends only on £ (and not m), an analogous 
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equation to (3) can be derived for cartesian harmonic coefficients: 

K<£l. ae (q) = 4tt| drr 2 K i (q,r)SS... ae (r) , (6) 

where the coefficients are defined as 

SSI ^r)J™±±» J ^A^U(^)S(r). (7) 
Cartesian harmonics satisfy the completeness relations: 

f (ff _ n) = ± E (H+M E ^ (!i0 ^ (n) 

= i E E >*2U ^ • • • • (8) 

^ i l - ai...at 

Replacing A^\VL) with the product n ai ■ ■ -n ae was justified by the fact that 
the projection operator is employed in A.^(£l'). Equation (8) allows one to 
restore 72. (q) and «S(r) from the cartesian coefficients, 

^(q) = E E ^L ae q ai ---q ae , 

5(r) = E E SS...« e r ai ---r ae , (9) 

£ ai...a e 

where q and r are unit vectors in the directions of q and r. 

Cartesian harmonics have other important properties. They are symmetric 
and traceless, i.e., if one sums «4^. over any two equated indices, the result 
is zero. In fact, the projection operator V in Eq. (5) is also referred to as a 
de-tracing operator. Tracelessness ensures, on its own, that r l A^ satisfies the 
Laplace equation. While the z axis is singled out in the construction of tesseral 
harmonics, the axes are treated symmetrically in the construction of cartesian 
harmonics, and different components can be obtained from each other by 
interchanging the x, y and z axes. For instance, the expression for A^ 3 ^ 
can be found by replacing replacing n x with n z in the expression for A^ 3 ^ 
in Table 1. Cartesian harmonics can be also generated recursively, beginning 
with ^=°) = 1, 
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•Aal...a e (ty ~ n X] ^ii..ii-ia, + i...«((^) 



£(2£ - 1) 



X ^ ^ ^ n Q a, +1 ...Q f (^) ■ (10) 



: i+ l...o;j_l Oj+i.-.a* v 

i<j a 



The above recursion can be used to obtain the coefficients for expanding the 
cartesian harmonics in terms of unit vector components [9]. The strategy of re- 
cursion, carried out in the other direction, leads to the coefficients of expansion 
of the products in terms of harmonics: 



m a <£ a /2 

_ (2€-4m + l)H „ £J 

Here, (—1)!! — 1, m — m x +m y + m z , £ — £ x + £ y + £ z and [4, 4, 4] is any set 
of indices where x occurs £ x times, y occurs £ y times and z occurs £ z times. 

The expansion (11) is useful for calculating moments of «S(r), 



f d 3 rx ex y iy z e * S(r) I [ d 3 r5(r) = V — — — 2m ^' — y**'*"'** 
J y u /i KJ /0 (2£-4:m + l)\\ Amx ' my > mz 

x / to^Sfc^_^_^lr) // drr 2 ^) . (12) 

In expansion in terms of cartesian harmonics, reflection symmetries with re- 
spect to the coordinate axes manifest themselves explicitly. Thus, if e.g. the 
source function is symmetric under the y — > — y transformation, all odd-£ y 
coefficients of S vanish. 

Since rank-£ cartesian harmonics are linear combinations of rank-£ tesseral 
harmonics, cartesian harmonics of different £ are orthogonal to one another. 
Within a given rank, though, there are {£ + 1)(£ + 2)/2 combinations of 
[£ x ,£y,£z], while only (2£ + 1) components can be independent. As an inde- 
pendent set, one can choose the components characterized by £ a = 0, 1, for 
a selected coordinate a, and the remaining harmonics within the rank can 
be generated via the traceleness property. For example, given harmonics for 
4 = 0, 1, the application of A[e x+2 ,e y ,e z ] = ~A^/v+^M ~ yields 
the other harmonics. As for scalar product values within a rank, they can be 
found in [9]. 

For weak physical anisotropies it is likely that only the few lowest terms of 
angular expansion can be ascribed to the correlation functions or the source. 
For £ < 2, the results can be summarized in terms of amplitudes and direction 
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Fig. 1. Characteristics of the sample source. The ellipse represents the contour of 
constant density within x-z plane, at half the maximum density. 

vectors as a function of q and r, as illustrated here with the case of TZ: 

ft(q) = K^\q) + J2 KS\q) q a + £ K$ a M *«i + - ■ (13) 



The function 7?.^ ' ) is the correlation averaged over ang les. The compo- 
nents describe the dipole anisotropy and can be summarized in terms of the 
direction e w (q) and magnitude R {1) (q), Tl£\q) = R {1) e£\ The TZ (2) compo- 
nents describe the quadrupole anisotropy and can be summarized in terms 
of three orthogonal principal direction vectors {e^\q)}k=i,2,3 and three dis- 
tortions that sum up to zero due to tracelessness: TZ^ai^) = 4 efi e ifl + 
e 3a e fp ~ + -^3 ) e 2a e i} '■ F° r higher-order distortions, the direction 
of maximal distortion elj may be defined as one maximizing the £'th carte- 
sian harmonic projected in the direction of e 3 : J2 ai ...a e ^-a}... ae e 3ai ' ' ' e 3a e ■ The 
condition of maximum implies E ai ...a e ft§ a 2 ...a e e kli 41 • • ■ e sL = 0, where 
k = 1,2 and ef and span the space of directions perpendicular to e^. 
The direction of may be chosen to minimize the projection of IZ^ in the 
plane perpendicular to e| . 

We next illustrate with an example how the correlation functions and source 
can be quantified following the cartesian coefficients. We choose a coordinate 
system with the z-axis along the beam direction and x-axis along the pair 
transverse momentum P -1 . We take a sample relative source for the particles, in 
a gaussian form with a larger 6 fm width in the beam direction and lower 4 fm 
widths in the transverse directions, cf. Fig. 1. The source is displaced from the 
origin by 5 fm, along P in system cm, assumed to point at 30° relative to the 
beam axis. The cartesian coefficients are calculated for the source according to 
Eq. (9) and expressed in terms of amplitudes and angles shown in Figs. 2 and 3. 
The displayed dipole characteristics result then from: sindg = S% = 
3/^<S(r)sin# and cos 6^ = Sjp = 3 / <S(r) cos and quadrupole 
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Fig. 2. Low-£ values of the sample relative source, as a function of the relative dis- 
tance r. Lines and symbols represent, respectively, exact values and values obtained 
through imaging of an anisotropic correlation function. Top panel shows the source 
averaged over angles, S^°\ Middle panel shows the dipole distortion . Bottom 

(2) (2) 

panel shows the larger, S% , and smaller, S| , quadrupole distortions within the 
x-z plane. 

from: S%> = S? ] cosHf + Sfsin 2 ^, Sg = sin 2 9® + Sf cos 2 ef 
and S?J = (fif - fif ) sin 2 where $£) = f J f 5(r) (sin 2 9 cos 2 - J), 

5® = f / f 5(r) (cos 2 - |) and Sg = f J f 5(r) cos (9 sin 5 cos 0. 

Despite the fact that the maximum of source density S occurs away from 
r = 0, the angle-averaged source in Fig. 2 is maximal at r = 0. An 
analytic source function S may be Taylor expanded in three dimensions about 
r = 0. The lowest-order terms of the Taylor expansion that can contribute to a 
rank-^ cartesian coefficient of S must involve an tth order derivative and rise 
as r £ , as is apparent in Fig. 2. The subsequent contributing terms rise as r e+2k , 
k = 1, 2, . . ., which is important if one tries to parameterize the functions S^ e \ 

At low r, the angles 9g and 9g are determined by the derivatives of S at 
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Fig. 3. Angles characterizing the dipole (top) and quadrupole (bottom) distortions 
of the sample relative source, as a function of the relative distance r. Lines and 
symbols represent, respectively, exact angles and angles obtained from imaging of 
an anisotropic correlation function. 

r = 0; in particular 6^ gives the direction of the gradient. At high r, the angles 
follow the elongation of the gaussian and, correspondingly, they approach zero 
as r — ► oo. 



Next, we determine correlation function coefficients for our sample source, 
examine how the coefficients reflect the source features and attempt to re- 
store those features through imaging. We choose the classical limit of repul- 
sive Coulomb interactions for the emitted particles, such as appropriate for 
intermediate mass fragments [12,13]. The kernel for Eqs. (1) and (4) can be 
analytically calculated by considering changes in momentum volume associ- 
ated with Coulomb trajectories, |0| 2 = d 3 q /d 3 q, where q is the starting 
momentum at the separation r. The result depends only on the angle # qr and 
on the separation at emission scaled with the distance of closest approach in 
a head-on collision, r/r c , where q 2 /2m a b = Z a Zi,e 2 /A7ieor c : 

|^)(g,r,cosMI 2 = 9(1 / +C ° Sgqr " 2 I c/r) ( 1 + COS ^- r -/ r ), (14) 

J (I + cos6> qr ) 2 - (1 + cos# qr ) 2r c /r 



The £ = kernel is K = 0(r — r c ) y 1 — r c /r — 1, while the Kt>\ components 
are, generally, calculated numerically from (4). 
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Fig. 4. A schematic relative source with repulsive Coulomb trajectories superim- 
posed coming out isotropically from position r within the source. The region of 
radius r c is inaccessible to the trajectories. The source consists of an isotropic 
part A, axially symmetric part Bi, responsible for source deformation, and, op- 
tionally, part B2 that is a mirror reflection of Bx- 

The Coulomb trajectories are illustrated in Fig. 4. The repulsive Coulomb force 
principally focuses the trajectories towards the direction of r they originate 
from. For ^ -C r (large q), the Kg kernels from from (4) are affected by the 
deflection of trajectories away from cos# qr ~ — 1. With Pi(cos6 — > — 1) = 
(— the kernels in this limit are Ki ~ (— l) e K Q ~ (— r c /2r, with 
the additional negative sign representing trajectory depletion. For r — r c <C 
r c (low q), the kernels are affected by the trajectories bunching up around 
cos# qr ~ 1. With Pi(cos8 — > 1) = +1, the kernels in this limit are K(>i ~ 

K + 1 = 0(r — r c ) \jr jr c — 1. The K\ kernel is always positive and, thus, the 
dipole distortion of the correlation generally points in the same direction as 
the distortion of the source, as for e.g. the source out of A and B x in Fig. 4. 
However, the K 2 kernel switches sign at r/r c = 1.63. At r/r c < 1.63, a prolate 
source, such as out of A and B 1 (and possibly B 2 ) in Fig. 4, results in a prolate 
correlation function elongated in the same direction. However, at r/r c > 1.63, 
the prolate source results in an oblate correlation with distortion pointing in 
transverse directions. 

Figure 5 shows the characteristics of the correlation function C, in terms of 
£ < 2 cartesian coefficients, for the source of Fig. 1, as function of r^ 1 / 2 oc q. 
The dependence 0^\r c ) in Fig. 5 generally retraces the dependence 0$~\r) 
of Fig. 3, where r c < r, as anticipated. However, while 9q (r c ) starts out in 
a similar fashion, it jumps by 90° as K 2 switches sign, with r in r/r c ~ 1.6 
characterizing the source quadrupole distortion. The jump is accompanied by 
the cross-over behavior for the distortions C^ 2 \ expected for the kernel sign 
change. In the low and high q limits, the correlation coefficients approach 
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Fig. 5. Low-£ characteristics of the Coulomb correlation function for our sample 
source: amplitudes (top panel) and angles (bottom panel) as a function of inverse 
square root of the head-on return radius r c . 



universal integrals of source coefficients, such as TZg) a (r c ) = Cf®^ a (r c ) — 
6io ~ 27rr c (— iy +l / °° drr a {r), for r c <C r. In that limit, from r^ 1 / 2 ~ 
0.9 fm™ 1 / 2 on in Fig. 5, the relative magnitudes of (— 1)^ reflect the relative 
magnitudes of the distortions for the integrated traceless tensors. 



We next illustrate imaging of source components from the correlation 
[10,14,15,16,17,18]. The tensorial decomposition allows to use different source 
representations within subspaces of different tensorial rank. With a source co- 
efficient represented in the basis {-Bj} as S(r) = J2i B^r), the 1Z function at 
j'th discrete value of q (or r~ 1 / 2 in our case), is IZj = KjiSi. Here, following 
(6), the matrix elements are Kji = 47T dr r 2 K(qj,r) Bi(r) and the spherical- 
tensor indices are suppressed. Minimization of x 2 — J2j(J^j ~ H ( f p ) 2 / v 2 with 
respect to {Si} yields the result, in matrix form, 



S = (K > ere 2 K)- L K T a c 2 TZ, 



(15) 
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where (c<J 2 )jfc = Sjk/cr 2 . Equation (15) is the basis of imaging, i.e. determining 
S^.. values from 7^^... Q ,. For illustration, we assume that the £ < 2 carte- 
sian coefficients of the correlation function for our sample source have been 
measured at 80 values of r^ 1 / 2 , between and 3fm~ 1//2 , subject to an r.m.s. 
error of a = 0.015. For the basis {Bi} we take take simple rectangle functions 
that, within basis splines, we generally find to provide the most faithful source 
values. As relative errors decrease with increasing multipolarity, we reduce the 
number of functions in restoration, from 5 for £ = to 3 for £ = 2, all span- 
ning the region r < 17 fm. Upon sampling the cartesian coefficients of 7Z, we 
restore the coefficients of S following (15) and obtain the results represented 
by symbols in Figs. 2 and 3 for central arguments of the basis functions. As 
is apparent, under realistic circumstances, angular features of the source can 
be quantitatively restored. 

Having the imaged £ < 2 source coefficients, we can calculate the £ < 2 
cartesian moments for the imaged region from (12). For the image we 
find, with statistical errors of restoration: 4nJdrr 2 S^ = 0.99 ± 0.05, 
(x) = (1/3) / drr 3 S^/ Jdrr 2 S { ^ = 2.47 ± 0.11 fm, (z) = 4.25 ± 0.13 fm, 
((x - {x)f) l/2 = 3.80 ± 0.24fm, {y 2 ) 1 ' 2 = 3.81 ± 0.22fm, ((z - {z)) 2 ) 1 ' 2 = 
5.54 ± 0.19 fm and ((x — (x))(z — (z))) = 2.23 ± 1.49 fm 2 , which can be com- 
pared to the r < 17 fm results from the original source of, respectively: 1.00, 
2.45 fm, 3.90 fm, 3.99 fm, 4.00 fm, 5.60 fm and -0.41 fm 2 . 

In summary, we have discussed the utility of cartesian surface-spherical har- 
monics in the analysis of particle correlations at low relative-velocities. The 
cartesian harmonics allow for a systematic quantification of anisotropic corre- 
lation functions, through expansion coefficients related to analogous expansion 
coefficients for anisotropic emission sources. For illustrating the relation, we 
have employed correlations produced by classical Coulomb interactions. To an 
extent, the features of source anisotropies may be read off directly from the 
correlation anisotropies; otherwise, they can be imaged. 
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